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ABSTRACT 

We numerically model the coronal mass ejection (CME) event of October 28, 2003 that erupted from 
active region 10486 and propagated to Earth in less than 20 hours causing severe geomagnetic storms. 
The magnetohydrodynamic (MHD) model is formulated by first arriving at a steady state corona and 
solar wind employing synoptic magnetograms. We initiate two CMEs from the same active region, one 
approximately a day earlier that preconditions the solar wind for the much faster CME on the 28th. This 
second CME travels through the corona at a rate of over 2500 km s^^ driving a strong forward shock. We 
clearly identify this shock in an image produced by the Large Angle Spectrometric Coronagraph (LASCO) 
C3, and reproduce the shock and its appearance in synthetic white light images from the simulation. We 
find excellent agreement with both the general morphology and the quantitative brightness of the model 
CME with LASCO observations. These results demonstrate that the CME shape is largely determined 
by its interaction with the ambient solar wind and may not be sensitive to the initiation process. We then 
show how the CME would appear as observed by wide-angle coronagraphs onboard the Solar Terrestrial 
Relations Observatory (STEREO) spacecraft. We find complex time evolution of the white-light images 
as a result of the way in which the density structures pass through the Thomson sphere. The simulation 
is performed with the Space Weather Modeling Framework (SWMF). 

Subject headings: MHD, coronal mass ejections, shock waves, numerical modeling 



1. INTRODUCTION 

The October-November period of 2003 saw some of the 
most energetic solar flares and coronal mass ejections of 
any solar cycle. These eruptive events came to be known 
as the Halloween Events and included 11 X-class flares, 6 
radiation storms and 4 geomagnetic storms. These erup- 
tions originated from three active regions (ARs), 10484, 
10486 and 10488 of which 10486 was the largest, most 
complex and most active producing 12 major coronal mass 
ejections (CMEs). Three of these CMEs that occurred on 
October 28, 29 and November 4 were particularly energetic 
with speeds in excess of 2000 km s~^ and were associated 
with X17, XIO and X28 flares respectively. The Octo- 
ber 28 CME was very geoeffective and damaged satellites, 
diverted airplane routes, caused power failures in Sweden, 
disrupted long-distance radio communications, and caused 
northern lights (aurora borealis) as far south as Florida. 
The November 4 CME was the most energetic of the Hal- 
loween Events, but occurred while AR 10486 was on the 
western limb so that the CME was less geoeffective as the 
CME was not directed toward Earth. 

Because of their extreme nature, the Halloween Events 
have been the subject of a great deal of study, both obser- 
vational as well as theoretical and modeling efforts. The 
October 28 CME in particular was well observed from 
Sun to Earth beginning wit h the vector magnetic field 
of AR 10486 (iLiu et al. l2005D. and surfac e flow measure- 
ments ()Yang et al.ll20o4lDeng et aI.ll2006D For this event , 
flare loops were observed by TRACE (|Su et all [20061 ) . 



and the global coronal field has been reconstructed from 
data taken from the Michelson Doppler Imager (MDI) 
aboard the Solar and Heliosphe ric Observations (SOHO) 
spacecraft ()Liu fc H avashi 2006V Dense plasma associated 
with the October 28 CME has been observed to propa- 
gate pas t the Earth by the Solar Mass Ejection Imager 
(SMEI) (iJackson et all [200l and was also detected b y 
interplanetary scintillation (IPS) (jTokumaru et al.ll2007[ ). 
Finally, the event was observed in situ b y Wind and 
the Advanced Co mposition Explorer (ACE) (jSkoug et al.l 
120041 : IZurbuchenI [2004), which reveal magnetohydrody- 
namic (MHD) quantities in the solar wind as well as ion 
abundances. 

Efforts to simulate the H alloween Events include those 
of llntriligator et al.l (|2005( ) who modeled the propaga- 
tion of the CME plasma from 2.5 Rs to 1 AU with 
a three-dimensional (3D) kinematic model. iKrall et"al] 
(|2006) modeled the 2003 October 28-30 period with a 
one-and-a-half-dimensional (1.5D) reduced MHD model 
that describes the self-similar expansion of a magnetic flux 
rope. In this case, the authors also coupled the helio- 
spheric results with a 3D magnetosphere code and com- 
pared itwith_smiu^^ by the observed solar wind 
data. iLiu &: Havashil (|2006f ) modeled the CME eruption in 
the solar corona with de nsity and p ressure perturbations. 
iDrver et al.l (|2004D and IWu et al.l (|2005D both modeled 
CME-driven shock propagation and predicted the shock 
arrival times with reasonable success. This work to model 
a specific event marks a clear distinction from non-event- 
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specific CME simulations (e.g. Usmano v fc Drvei 19951: 
Rilev et all 120021 : [Manchester et al. 2004bt [Jacobs et all 



ey I 

2003), ELnd non-event-specific Sun-to-thermosphere space 



weather simulations modeled b y the Center for S pace En- 
vironment Modehng (CSEM) i|T6th et al.l 120051) and the 
Center for Integ rated Space Weather Modeling (CISM) 
(jLuhmann e t al.l l2004D. 

More recentlv. lToth et al.l (|2007D modeled the October 
28 CME as part of a Sun-to-Earth space weather event 
that included an in-depth des cription of the resu lting geo- 
magnetic storm. This work bv lToth et al.l (|2007) is ground 
break i ng in two regards. F irst, hke those by iLugaz et al.l 
(|2007t ): ICohen et al.l (|2007t) . it is the first 3D numerical fuU 
MHD simulation to reproduce observed CMEs as they are 
magnetically driven from active regions in the low corona. 
The earliest attempt to model an observed CME (includ- 
ing the magneti c field) propagat ing from the Sun to the 
Earth, is that of lWu etal] (|l999l ). This particular simula- 
tion is two-dimensional and did not capture the structure 
of the solar wind or the active region, but did roughly 
capture the magnetic cloud at 1 AU. More recent mod- 
eling efforts have been 3D and capture the structure of 
the solar wind, but treat CME propagation outside of the 
magnetosonic point (r « 20Rs) (jOdstrcil et al.l 20051). The 



second noteworthy aspect of the simulation bv lToth et al.l 
(|2007( ) is that it is the first event study to employ a frame- 
work to couple several individual codes to model the physi- 
cal domain extending from the corona to the Earth's upper 
atmosphere. The framework used was the Space Weather 
Modeling Framework (SWMF) develo ped by members o f 
CSEM at the University of Michigan bv lToth et all (|2005D . 
A framework for similar purposes is being developed by 
CISM at Boston University. 

In this p a per, w e examine the simulation described by 
iToth et al.l (|2007D . and concern ourselves with the phys- 
ical properties of the CME as it propagates from the 
low corona to Earth. In particular, the detailed struc- 
ture of the coronal and solar wind allow us to quanti- 
tatively compare the model with Large Angle Spectro- 
metric Coronagraph (LASCO) C2 and C3 observations of 
the October 28 event. In this case, we validate the ac- 
curacy of the simulation in reproducing the speed mass, 
appearance (in scattered light), and shock properties of 
the observed CME. We also make wide-angle, large elon- 
gation Thomson-Scattered white light images of the model 
to show how the CME would appear in the wide angle 
coronagraphs of the Sun Earth Connection Corona and 
Heliospheric Imager (SECCHI) instrument on board the 
STEREO spacecraft. These synthetic images show that 
the CME appears much different far from the Sun due to 
scattering affects at large elongation. 

The rest of this paper is organized as follows: The 
SWMF, its components, and numerical techniques are 
briefly discussed in section 2. Section 3 provides a de- 
scription of the CME simulation, the results of which are 
compared to LASCO observations in section 4. In sec- 
tion 5, we show how the CME would appear propagating 
past the Earth in the SECCHI heliospheric imagers HIl 
and III2. We conclude with a summary and an outlook 
for future development in section 6. 

2. NUMERICAL METHODS: SWMF AND BATS-R-US 



The SWMF (|T6th et al.l I2005D couples models treat- 
ing physical domains of the space environment extending 
from the solar corona to the Earth's upper atmosphere. 
The model coupling is flexible yet efficient, making faster 
than real-time space weather simulations feasible on mas- 
sively parallel computers. Each model has its own de- 
pendent variables, a mathematical model, and a numer- 
ical scheme with an appropriate grid structure and tem- 
poral discretization. The physics domains may overlap 
with each other or they can interact through a bound- 
ary surface. The SWMF is able to incorporate models 
from the community and couple them with modest changes 
in the software of an individual model. The SWMF is 
a fully functional and documented framework that pro- 
vides high-performance computational capability to sim- 
ulate the physics from the low solar corona to the upper 
atmosphere of the Earth. Currently, the SWMF is com- 
posed of nine physics modules, and is driven by external 
data such as magnetograms, flare and CME observations, 
satellites upstream of the Earth (like ACE, Geotail and 
Wind), etc. The SWMF can model any physically mean- 
ingful subset of the physics domains, and is freely available 
via registration at 

http : / / csem . engin . umich . edu/SWMF' 

There were seven SWMF components used in the Hal- 
loween storm simulations: Solar Corona (SC), Eruptive 
Event Generator (EE), Inner Heliosphere (IH), Global 
Magnetosphere (CM), Inner Magnetosphere (IM), Iono- 
sphere Electrodynamics (IE), and Upper Atmosphere 
(UA). For our purposes, the only relevant components are 
the SC, EE and IH that model the corona, solar wind, and 
CME propagation from the active region to 1 AU. Com- 
pone nts SC, EE, and IH are based on the BATS-R-US 
code (jPowell et al.lll999l: iGombosi et al.ll200lD that solves 
the governing equations of MHD in a conservative finite 
volume form. Non-ideal MHD terms are included through 
appropriate source terms. The code uses a limited recon- 
struction that ensures second-order accuracy away from 
discontinuities, while simultaneously providing the stabil- 
ity that ensures non-oscillatory solutions. In addition, the 
code employs several approximate Riemann solvers. The 
resulting scheme solves for the hydrodynamic and electro- 
magnetic effects in a tightly coupled manner that works 
equally well across several orders of magnitude in plasma 
P (the ratio of plasma pressure to magnetic pressure). 
The BATS-R-US also uses a relatively simple yet effec- 
tive block-based adaptive mesh refinement (AMR) scheme 
to resolve structures spanning many length scales. The 
BATS-R-US scales almost linearly to more than 1000 pro- 
cessors for a fixed problem size, and fully implicit time 
stepping scheme is incorporated, (|T6th et al.ll2006l ) that 
can be combined with explicit time stepping. 

The Solar Corona domain is a Cartesian box that ex- 
tends from the surface of the Sun to — 24Rs < x,y,z < 
24 Rs , where Rs is the radius of the Sun. The physics of 
this domain is described by the equations of MHD with 
additional source terms required to take into acc ount the 
heating and acceleration o f the solar wind (Groth et aLl 
120001: lUsmanov et al.|[20M) . Recently. ICohen et al.l (|2007^ 
developed a solar wind model based on the empirical rela- 
tionship be tween solar wind speed and magnetic flux tube 
expansion (|Wang fc Sheelevlll990l) . which has been mcor- 
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Fig. 1. — The initial condition of the corona for the October, 28 CME. The panel on the left (a) shows the radial field strength, Br at the 
base of the corona. The structure of the coronal magnetic field is illustrated with blue and yellow lines for the closed field while orange lines 
show the open field. The panel on the right (b) shows a close up of AR 10486 with the superimposed flux rope illustrated with field lines and 
a current density iso-surface colored green. 



porated into the SWMF. Here, we use the coronal model 
presented in iRoussev et al.l (|20031j ) . At the inner bound- 
ary of the SC component, the density, pressure, velocity 
and magnetic field arc defined at a height just above the 
transition region. The magnetic field is obtained from a 
synoptic solar magnetogram. The boundary conditions for 
the temperature and mass density at the Sun are varied as 
a function of magnetic field strength to achieve a realistic 
distribution of fast and slow wind speeds near the Sun and 
at lAU. The velocity components at the inner boundary 
maintain line-tying of the magnetic field to the rotating 
solar surface. Differential rotation is currently neglected. 
The flow at the outer boundary is usually superfast (faster 
than the fast magnetosonic speed of the plasma), so no in- 
formation is propagating inward. 

The EE domain is in the Solar Corona, restricted to the 
active region responsible for the CME. The EE in this case 
takes the form of a nonlinear perturbation of the SC solu- 
tion, which is made by superirn posing a modified version 
of the lTitov fc DemoulinI (|l999t) flux rope to active region 
10486. The flux rope is in an initial state of force im- 
balance that drives the eruption. This eruption generator 



follows from t he work of iRoussev et al.| (j2003a[ ) who in- 

30) flux rope in the 



corporated the iTitov fc Demoulin ( 199E 
BATS-R-US code. 

The IH domain extends from its inner boundary at r = 
20 Rs to a cube extending to -240 Rs < x,y,z < 240 Rs, 
which encompasses Earth's orbit. The physics of this do- 
main is described by the equations of ideal MHD, solved in 
an inertial frame, on a Cartesian grid. The inner bound- 
ary conditions of the IH component are obtained from the 
SC component. The flow at the outer boundary of the IH 
component is always assumed to be superfast. 

3. HALLOWEEN STORM SIMULATION 

3.1. Initial State 

Here, we summarize the basic desc ription of the CME 
simulation that was first presented in iToth et al.l (|2007[ ). 



For more details, we refer the reader to this earlier paper. 
The simulation begins with the construction of a steady 
state corona and solar wind. The magnetic field at the 
inner boundary (base of the corona) of the SC component 
is specified from an MDI synoptic magnetogram centered 
around the time of the October 28 eruption. The initial 
volumetric field is specified with a potential-field-source- 
surface extrapolation fit to this map. The computational 
grid is highly refined around AR 10486 where the smallest 
cells are about 3 x lO^'^Rs. The grid is refined to 0.1 Rs 
within a 1 Rs diameter cylinder extending along the Sun- 
Earth line. The SC component is then run to steady state 
allowing the solar wind to relax. Figure[l] Panel (a) shows 
the steady state solution. Here the radial field strength 
is shown in color at the inner boundary, and field lines 
are drawn extending out into the corona. Orange lines 
are open while closed lines are shown in blue. The SC is 
run in the Heliographic Rotating (HGR) coordinate sys- 
tem with the sidereal Carrington rotation period of 25.38 
days. An offset angle around the Z axis places the Earth 
in the —X, Z half plane at the time of the large eruption 
on October 28 CME. 

The IH component is coupled to the SC component at 
20 Rs, and rapidly achieves steady state as the superfast 
wind blows through the domain. During the steady state 
run, the IH component uses a rotating coordinate system, 
which is then switched to Heliographic Inertial (HGI) with 
an offset angle that puts the Earth in the —X, Z half plane 
(the orbital motion of the Earth is neglected). The grid 
resolution is 1/4 Rs near the inner boundary and also along 
the Sun-Earth line within a cylinder of 1 Rs radius, while 
the largest cells are 4 Rs . 

3.2. Generating Eruptive Events 

The October 28 event was preceded by several smaller 
CMEs, which significantly modified the ambient solar 
wind. To take into account the preconditioning effect of 
previous CMEs on the solar wind, we start our Simula- 
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Fig. 2. — The structure of the CME at time after initiation (TAI) equal to 15 minutes and 45 minutes is shown in the top and bottom rows 
respectively. In Panels (a) and (c) the electron density is shown in color with magnetic stream lines (confined to the plane) drawn white and 
the numerical mesh drawn black. These images show the ejected magnetic flux rope traveling toward the Earth in the —x direction with the 
center of the rope and densest plasma concentrations below the equatorial plane. Panels (b) and (d) show the ratio of the density relative 
to the pre-event state on j; — 2 planes positioned near the center of the CME at x = — 5Rs and x = — 8Rs, respectively. The white and 
black circles are drawn at 2 and 3.7 Rs, respectively, corresponding to the occulting disks of C2 and C3, respectively. These images show the 
greatest density enhancement at the bottom (— z) of a nearly circular shell type structure that evolves in a self-similar fashion. 



tion on October 26, when a smaller CME was launched 
at approximately 07:00 UT from AR 10486. For this first 
CME, we use 07:24 UT as the initiation time, which is 
within the uncertainties of the observations. The much 
more energetic CME that is the focus of this study oc- 
curred at approximately 10:54 UT. This eruption occurred 
near disk center producing a full halo CME. 

Here, we initiate the CMEs by inserting magnetic flux 
rop es (the size of the a ctive region) based on the modi- 
fied iTitov fc DemoulinI (|1999D model. In Figure [1] (Panel 
(b)), the flux rope for the October 28 eruption is seen 
arching above AR 10486. The location and orientation of 



the flux ropes were chosen to arch over the dipolar part 
of the active region. The active region can be easily iden- 
tified in the high resolution synoptic map obtained from 
the observed photospheric magnetic field. The density in 
the loop was obtained from the size and estimated mass 
of the CMEs. The rope's magnetic field does not match 
that observed at the photosphere. The reason for the dis- 
crepancy is that the field of the rope is highly twisted 
and is orthogonal to the polarity inversion line. In con- 
trast, the measured magnetic field was observed to be in a 
highly sh eared configurat ion running parallel to the inver- 
sion line ()Liu et al.l 2005f ). In spite of this limitation, the 



MANCHESTER ET AL. 



5 




Fig. 3. — Three-dimensional density structure of the model CME 15 minutes after initiation. Iso-surfaces of the density enhancement 
(relative to background) are shown at levels of 0.9, 2 and 5 in blue, green and red respectively. The green iso-surface shows the extent of the 
shock front, while the blue iso-surface shows the low density rarefaction behind the CME. The Sun is shown with a yellow sphere. 
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Fig. 4. — Comparison of observed and modeled CME velocity. The solid line shows the modeled CME velocity moving directly toward the 
Earth, while the dashed line shows the model velocity projected on the plane of the sky 177 degrees (counter clock wise) from the north polar 
axis. At this same location in the plane of the sky, the CME velocity is derived from LASCO observations and plotted with stars. We find 
that the model briefly reaches a velocity of 4000 km s~^ at 4.5Rs before falling to 2000 km s~^ at 20Rs. In contrast, the CME is observed 
to decelerate from 1890 km s~^ to 1699 km s~^ as it travels from 2.3Rs to 26.3Rs. 



magnetic field strength of the flux ropes can be adjusted to produce eruptions that match LASCO observations of 

the CME speeds at 20 solar radii (1500 and 2500 km s^^ 
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Fig. 5. — Comparison of observed (left) and simulated (right) Thomson-scattered white light brightness. The color images show the total 
brightness divided by that of the pre-event background. White circles show the solar limb, filled black circles show occulting disks. Panel (a) 
shows the LASCO C2 observation at 11:30 UT, while LASCO C3 observations are shown in panels (c) and (e) at 11:42 and 12:18, respectively. 
Panels (b) (d) and (f) show model coronagraph images at 13, 21 and 50 minutes after initiation. Here, we find that the model achieves 
remarkable quantitative agreement with both the magnitude and spatial distribution of the observed brightness. 



for the October 26 and October 28 events, respectively). 
For the first CME, the magnetic flux rope has a free en- 
egy 2.3 x 10^^ ergs, while the second faster CME requires 
nearly ten times as much energy at 2.0 x fO'^^ ergs. 

After 20 minutes from the start of the time dependent 
simulation, the leading shock of the first CME reaches 
a radial distance of 5.5Rs at a speed slightly exceeding 
2100km s~i. The first CME reaches the SC/IH bound- 
ary at 20 Rs after about 1.7 hours, and the speed of the 
leading shock is the observed 1500 km s^^. The first CME 



is propagating in a direction about 30 degrees off from 
the Sun-Earth line, but the flanks of the shock reach the 
Earth 45 hours after the eruption at 4:30 UT Oct 28. The 
simulated solar wind velocity at the Earth increases from 
about 350 km s~^ to 550 km s~^, which is a good approx- 
imation to the solar wind conditions preceding the arrival 
of the October 28 CME as observed by the ACE satellite 
(|Skoug et al.ll2004l : IZurbuchenll2004D . 

The second CME is initiated at 10:54 October 28 with 
a flux rope that is 50% larger in radius than the previous 
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Fig. 6. — Shock structure as revealed by coronagraph images. Panels (a) and (b) respectively show color images of the total brightness 
(divided by the pre-event brightness) for LASCO C3 (time = 12:18) and for the simulation 50 minutes after initiation. The color legend is 
chosen to highlight the faint (2 %) increase in brightness at the shock found at the outer edge of the halo. The model shows the same faint 
increase in brightness at the shock as well as the capturing the width of the faint plasma sheath. Panels (c) and (d) show respectively line 
plots of the brightness and electron density as functions of distance from Sun center along the black line at 177 degrees from the north polar 
axis. Model and observations are shown with dashed and solid lines respectively. The line plots reveal that the simulation quantitatively 
reproduces the observed brightness increase behind the shock as well as the inferred density increase at the shock. 



one. Due to the larger size and stronger magnetic field, 
the second CME reaches a speed of 2500 km s~^ (on the 
Sun-Earth line) measured at 20 Rs. After 15 minutes, the 
second CME reaches 9 Rs with a shock speed near 3200 km 
s~^. The shock reaches the SC/IH boundary at 20 Rs in 
less than an hour with a speed around 2800 km s~^. The 
density structure of the October 28 CME is shown in Fig- 
ure m Panels (a) and (c) show a color representation of 
the electron density on the meridional plane {x — z) at 15 
and 45 minutes after CME initiation, respectively. White 
lines show the direction of the magnetic field confined to 
the plane. In these pictures, the magnetic flux rope is 
clearly seen to be expelled from the corona. Panels (b) 
and (d) show the density ratio (relative to the pre-event 
corona) on y — z located at x = — 5Rs and x = — BRs 
at 15 and 45 minutes after initiation, respectively. The 
white and black circles correspond to the solar limb and 
LASCO C3 coronagraph, respectively. These images show 
the nearly self-similar evolution of the plasma expelled in 
the CME. The plasma is distributed in a shell with the 
greatest density in the southern hemisphere. 

The three-dimensional density of the CME is shown in 
Figure[3]with iso-surfaces of the density enhancement (i.e. 
density divided by the pre-event level) at t = 15 minutes 



after initiation. Blue, green and red iso-surfaces corre- 
sponding to values of 0.9, 2, and 5 are shown from four 
perspectives: (a) side, (b) polar, (c) back and (d) front fac- 
ing Earth. The green surface effectively shows the extent 
of the shock front, the red surface shows the core mate- 
rial, while the blue iso-surface shows the density depleted 
cavity. The cavity forms as a rarefaction behind the flux 
rope, while the densest plasma is located in the southern 
hemisphere, both in and below the flux rope. The polar 
and back views (Panels (b) and (d) respectively) best show 
the density enhancement at the southern en d of the CME. 
Comp arison of Figure O with Figure 6 in iJackson et alj 
(|2006l ) also suggests that the simulated 3D density struc- 
ture qualitatively agrees with the density enhancements 
reconstructed from SMEI observations. 

4. COMPARISONS WITH LASCO OBSERVATIONS 

CMEs are most frequently observed in visible light that 
is Thomson scattered by coronal electrons within the 
CME. In addition to this electron scattered photospheric 
light (referred to as the K-corona), there is also light scat- 
tered by interplanetary dust that is referred to as the F- 
corona. There are significant differences between the two 
components of scattered light. First, the F-corona is unpo- 
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larized within 5-6 Rs, while the K-corona is polarized with 
components b oth radial and tangential to the solar limb 
("Billings! [19661) . Second, the brightness of the K -corona 
falls off much more rapidly with distance from the Sun 
such that it dominates close to the Sun, is roughly equal 
to the F-corona at r ss 4Rs, and is much dimmer than the 
F-corona far from the Sun. Finally, the F-corona is nearly 
constant in time and is believed to be largely unaffected 
by CMEs. 

The corona is optically thin in white light, so that the 
intensity of scattered light along a given line of sight (LOS) 
is the integrated contribution of all electrons along the line. 
The contribution to a given LOS is highly dependent on 
the locati on of the elec tron relative to the Sun and the 
observer (|Biningslll966l ). Close to the Sun, the scattering 
is heavily weighted in the plane of the sky. Coronagraphs 
located in space observe coronal light without contending 
with atmospheric scattering, which makes it possible to 
view the corona far from the Sun. In the case of LASCO, 
the C2 and C3 coronagraphs have fields of view extend- 
ing to 6 and 32 Rs respectively. More recently, SMEI and 
STEREO coronagraphs have observed Thomson scattered 
light at very large angles from the Sun that extend beyond 
Earth's orbit. 

In this section, we make comparisons between synthetic 
coronagraph images constructed from our model CME, 
and LASCO observations. The synthetic images are cre- 
ated by numerically integrating the Thomson scattered 
light along a LOS for each pixel with the appropriate scat- 
tering function (Billings 1966). Images of nume rical CME 



models have been made in this fashion before (IWu et al 



19991: IOdstrcini2003HManchester et alll2004at ITugaz et al 



20051 e.g.). More recently, synthetic images have been 
made for models of specific CME ev ents, and qualita- 
tively compared to LASCO images (e.g. lLugaz et al.ir2007t 
ICohen et al.ll2007fl . Here, we go further, and for the first 
time make truly quantitative comparisons with the data. 
To make the comparison, we process the synthetic data 
and the LASCO data in the same way and display the 
resulting images in identical formats. We account for F- 
corona far from the Sun where it is the dominant source of 
scattered light. Since our MHD model only allows us to di- 
rectly treat electron density, we estimate the contribution 
of dust scattering with the same power-law empirical re- 
lationship used to make the LASCO images. In this case, 
the background F-corona brig htness is taken to have th e 
form Bp = cr(°-22 ™«(2e)-2.47) (|Koutchmv fc LamvlfigSSh . 
Here, we estimate the magnitude of c by setting Bp = Bk 
at r = 4 Rs where Bk is th e background pola r brightness 
of the K-corona taken from ISaito et"al] ()1977[ ) . 

The results of the comparison are shown in Figure [5] 
Here, images of the October 28 CME as observed with 
LASCO C2 (Panel (a)) and C3 (Panels (c) and (e)) are 
shown for times t = 11:30, 11:42 and 12:18 respectively. 
Panels (b), (d), and (f) of Figure [5] show the synthetic 
images at times t — 13, 21, and 50 minutes after initia- 
tion when the model CME is closest in size to that of the 
corresponding LASCO image. In reality, the CME travels 
slower than in our model with the observed sequence of 
images occurring at approximately 36, 48, and 84 minutes 
after initiation. In both data and model, the color images 
show the total brightness divided by that of the pre-event 



background. White circles show the solar limb, and filled 
black circles show the occulting disks. In all cases, we find 
extremely good quantitative agreement in both the mag- 
nitude and spatial distribution of the observed brightness. 
In these images, the greatest brightness is below the south- 
ern limb of the Sun, that extends in an arc over the eastern 
limb, while the the dimmest part of the CME is found in 
the north western quadrant. This brightness distribution 
corresponds with the position of AR 10486, which was lo- 
cated at 20 degrees below the equator on October 28. By 
12:18 UT, the observed brightness patterns is more highly 
structured than the model, showing evidence of a bright 
core in the south. The model (Panel (f)) lacks these fine 
details, and retains a single bright arc in the south. 

4.1. CME Mass 

To further quantify the comparison between the our nu- 
merical model and the LASCO coronagraph observations, 
we calculate the CME mass in both cases in an identical 
fashion. This mass is derived by integrating the excess 
brightness of the K-corona with the assumption that the 
plasma is in the plane of the sky, which sets a lower limit 
on CME mass. The F-corona is subtracted out from the 
observations and does not affect the mass. We restrict this 
integral to the brightest portion of the CME that extends 
from 235 to 330 degrees as measured from the —y axis 
in Figure [5l In the case of the C2 field of view (top row 
of Figure [5]), this integration sector extends from 2.0 to 
6.0 Rs, and for the C3 field of view (middle and bottom 
rows of Figure [5]) the sector extends from 4.0 to 16.4 Rs. 
With these criteria, the CME masses derived from obser- 
vations at t = 11:30, 11:42 and 12:18 UT are 1.50 x lO^*^ g, 
1.55 X 10^® g, and 1.70 x 10^^ g, respectively. The corre- 
sponding masses of the model CME (Panels (b), (d), and 
(f)) are, respectively, 2.23 x 10^^ g, 2.16 x 10^*^ g, and 
4.90 X 10^6 g. 

The observed masses are approximately 30% less than 
that of the model at times t = 11:30 UT and 11:42 UT, 
which is consistent with the larger filling factor of the 
model event. The slight decrease in model CME mass 
found from Panel (d) is the result of a greater part of 
the CME being obscured by the occulting disk. By 12:18 
the CME mass is observed to increase by 13 %, while the 
model mass (as shown in Panel (f)) more than doubles. 
We expect the mass of fast CMEs to increase as plasma 
is swept up as they travel from the Sun as shown in nu- 
merical simulations (jManchester et al.ll2004al : iLugaz et al] 
I2005f) . However, the mass increase seen in LASCO obser- 
vations is attributed to mass ent ering the coronagraph field 
of view from below the occulter (IVourli das et al.l 20001 ). In 
general, there is not yet proof of a plasma pileup in the 
low corona by CMEs. In the case of our simulation, we 
believe that the CME mass increases faster than what is 
observed because of the excess speed of the CME. For the 
observed event, the (excess) mass continues to increase by 
almost an order of magnitude, to 13.6 x 10^^ g as measured 
with SMEI on October 29, 12: 00 UT as the CME passed 
the Earth () Jackson et al.ll2006f) . 

4.2. Shock Identification 

Faint arcs are frequently observed at the outer edges 
of fast CMEs viewed in coronagraph images, yet it has 
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Fig. 7. — Views of the model CME as it would appear in the SECCHI coronagraphs. The CME is shown in a time series propagating 
through the fields of view of CORl, COR2, HIl and HI2 that extend from 1.4 to 4Rs, 2 to 15Rs, 12 to 84 Rs (20 degrees wide), and 66 to 
381 Rg (70 degrees wide) respectively. In this case, the HI images are taken from the point of view of STEREO A, 22 degrees ahead of the 
Earth, while COR images are made from the location of Earth. 



been difficult to conclusively identify these arcs as shocks. 
The earliest exampl e of s hock identification is that of 
ISime fc HundhausenI (jl987f ) who observed a bright loop 
at the front of a fast CME. In that case, the presence of 
a shock was inferred from its high speed (1070 km s^^), 
the absence of any deflections preceding the shock, and 
the fact that the expanding shock front did not cease its 
lateral motion t o form stationary legs. More recently, 
IVourlidas et al.l (|2003f ) observed a similar feature with 
LASCO that they identified as a shock wave. In this work, 
a numerical simulation was employed to model the shock's 
appearance and confirm its presence in the coronagraph 
image. More recent MHD simulations have also found 
manifestations of C ME driven shocks in synthetic corona- 
graph images (e.g. [Manchester et al.l l2004at iLugaz et al.l 
|2007() . Perhaps the most compelling observational evi- 
dence for shocks appearing as visib le components of CMEs 
in L ASCO images i s prese nted in iRavmond et al] (|200Cl[ ) 
and iMancuso et al.1 (|2002[ ). Here, shocks were observed 
simultaneously in the low corona (r < 3Rs) by LASCO, 
the Ultraviolet Coronagraph Spectrometer (UVCS) and 
as type II radio burst. UVCS gave clear spectroscopic ev- 
idence for the presence of shock fronts, while radio bursts 
indicated the presence of shock-accelerated electrons. 



In our simulation, we can conclusively identify a shock 
front and determine its appearance in synthetic corona- 
graph images. By detailed comparison with the observa- 
tions, we can then clearly demonstrate the existence of a 
shock front in the corona. For this purpose, we reexamine 
the coronagraph data presented in the bottom row of Fig- 
ure [5l We again show the brightness, but with the color 
contours adjusted in Figure [6] to more clearly show the 
faint features at the outer-most edge of the observed and 
modeled CME in Panels (a) and (b) respectively. In both 
panels, we see very similar faint rims at the outer edge of 
the CME where the color makes a transition from blue- 
green and yellow indicating an increase in brightness at a 
level between 1 and 2%. This rim has a width of 2 Rs and 
extends out to 16 Rs from the Sun. In the data, this faint 
rim is near the noise level and we draw a curved black line 
to mark its outer boundary. 

For a more detailed comparison, we extract the bright- 
ness and electron density along radial lines (shown in Fig- 
ure [H Panels (a) and (b)) that pass through the southern 
outer edge of the CME. Panel (c) shows both the observed 
and modeled brightness with solid and dashed lines respec- 
tively. We find very good agreement between the model 
and the observations. Both cases exhibit a similar steep 
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Fig. 8. — Synthetic image of the CME as it would appear from SECCHI-HI2A one year into the mission (22 degrees ahead of the Earth). 
The total brightness is shown in a time sequence at intervals t =, 9.1, 13.1, 18.3, and 22.3 hours in panels (a), (b), (c) and (d), respectively. 
Here, we see the CME propagate into the field of view on the right hand side of Panel (a), appear to stall in Panel (b), and then begin to 
brighten on the left hand sides of Panels (c) and (d). This type of evolution is not observed close to the Sun, and occurs at large elongation 
because of the passage of the CME density structure over the spherical surface of maximum scattering. 



drop in brightness that levels off forming the faint rim, 
where the noise in the observations is close to the sig- 
nal level. For the purpose of the comparison, we have 
aligned these features in the radial direction by shifting 
model results a small amount (+0.62 Rs). In Panel (d), 
we plot (solid line) the electron density deduced from ob- 
servations under the assumption of an axisymmetric dis- 
tribution around the Sun. The model electron density is 
extracted along a line extending from the center of the 
Sun, and passing through the southern most extent of the 
disturbance near x,y,z = —8.6, 0.83, —15.9 Rg. The pro- 
jection of this line on the y — z plane corresponds to line 
shown in Panel (b). We plot this data in Panel (d) with a 
dashed line and translate the model density in the radial 
direction to align with the observed quantities. We find 
remarkable agreement between the two electron densities, 
where both increase from approximately 2000 to 10000 



cm^^ at the outer edge of the faint rim between 14 and 
ITRg. In the simulation, we can conclusively identify this 
density jump as a strong fast-mode shock driven by plasma 
moving at 2010 km s~^. The density increase at the shock 
is very near the theoretical limit of (7 -|- l)/(7 — 1) = 5 for 
7 = 1.5 used in the simulation. The shock is smeared out 
over a distance of 2Rg by grid resolution at 1/4 Rs. By 
a detailed comparison between the simulation and the ob- 
servations, we clearly demonstrate the existence of a shock 
front in the LASCO C3 image at t = 12:18. 

The overall 3D shape of the shock front is seen in Figure 
[3l from which we can discern why the shocked plasma is 
so faint. The outer edge of the shock front is 5Rs from 
the plane of the sky and is sharply curved at its greatest 
distance from the x axis. The lines-of-sight at the outer 
edge of the CME graze the shock surface, and even though 
there is nearly a factor of four increase in density at shock. 
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Fig. 9. — Synthetic image of the CME as it would appear from SECCHI-HI2B one year into the mission (22 degrees behind the Earth). 
The total brightness is shown in a time sequence at intervals t =, 9.1, 13.1, 18.3, and 22.3 hours in panels (a), (b), (c) and (d), respectively. 
Here, we see the CME propagate into the field of view on the left hand side of Panel (a), noticably brighten in Panel (b), fade in Panel (c), 
and then all but disappear in Panel (d). This evolution is much different than that shown in Figure[8] and illustrates the 3D structure viewed 
in the opposite direction as seen by the HI2A coronagraph. 



the enhancement occurs over such a short distance that it 
resuhs in a very small increase in brightness. The bright- 
ness increases as the lines of sight pass through progres- 
sively longer distances of compressed plasma behind the 
shock. For the numerical model, the decrease in bright- 
ness occurs faster than what is observed, in part because 
the shock front is smeared out by the limited numerical 
resolution. 

5. CME APPEARANCE AT LARGE ELONGATION: 
STEREO-SECCHI PREDICTIONS 

In this section, we make wide-angle, large elongation 
Thomson-scattered white light images of the model to 
show how the CME would appear in the coronagraphs of 



the Sun Earth Connection Corona and Heliospheric Im- 
ager (SECCHI) inst rumen ts onboard the STEREO space- 
craft (Ho ward et all 120071 ). The fields of view of the four 
SECCHI coronagraphs: CORl, C0R2, HIl, and HI2 give 
complete coverage extending from the low corona to be- 
yond the Earth's orbit as shown in FigureEl CORl, C0R2 
coronagraphs are centered on the Sun with fields of view 
that extend from 1.4 to 4Rs and 2 to 15 Rs respectively. 
HIl and HI2 are pointed toward the Sun-to-Earth line 
with fields of view that extend from 12 to 84 Rs and 66 
to 381 Rs, respectively. Figure [7] shows synthetic CORl 
and C0R2 images of our model of the October 28 CME 
as it appears from the location of the Earth. HIl and HI2 
images show the CME from a location 22 degrees ahead 
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Fig. 10. — The time evolution of the density structure of the CME on the equatorial plane. The density is shown in color in a time sequence 
at intervals t =, 9.1, 13.1, 18.3, and 22.3 hours in panels (a), (b), (c) and (d) respectively. The locations of the Sun, Earth, and STEREO A 
and it's field of view are shown along with the equator of the Thomson sphere, the location of maximum scattering into the observer's LOS. 
The panels correspond to the same times shown in Figure [S] and clearly show how dense plasma crosses the Thomson sphere to produce the 
time evolution of brightness seen in Figure \8\ 



of the Earth where STEREO A will be one year into it's 
mission. 

Thomson scattering has a strong angular dependence 
such that scattering into a LOS is highly dependent the 
on loc ation of the ele ctron relative to the Sun and the ob- 
server (lBillingslll966h . The region of space for which there 
is maximum scattering into LOS is in the shape of a sphere 
that contains the Sun and the observer as see in Figure 
[TUl Close to the Sun, the scattering is heavily weighted 
to the plane of the sky. At large elongation (more then 
60 Rs from the Sun), the spherical shape of the maximum 
scattering surface beco mes significant to the app earance 
of coronagraph images (|Vourlidas fc Howardl[200i) . 

Our model provides an opportunity to examine the ef- 
fects of the angular dependence of Thomson scattering on 
the appearance of CMEs at large elongation as seen in the 
HI2 coronagraphs. Figure [5] shows the time evolution of 
the model CME as it would appear from SECCHI-HI2A 
one year into the mission (22 degrees ahead of the Earth). 



The total brightness is shown in a time sequence at inter- 
vals t 9.1, 13.1, 18.3, and 22.3 hours in panels (a), (b), 
(c) and (d), respectively. The CME is seen to propagate 
into the field of view on the right hand side of Panel (a) , 
appears to stall in Panel (b), and then reappears on the 
left hand sides of Panels (c) and (d) without ever being 
bright in the middle of the field of view. (A movie is avail- 
able for the online version of this paper.) Figure [9] shows 
the corresponding series of images that would be made be 
SECCHI-HI2B. In this case, the spacecraft is 22 degrees 
behind the Earth with the HI2B instrument pointing to- 
ward the Sun-Earth line in the opposite direction as HI2A. 
From this perspective, the CME enters the HI2B field of 
view from the left hand side seen in Figure [9l Panel (a). 
Here, the CME appears as an arc that first brightens as 
seen in Panel (b) and then gradually fades away as seen 
in Panels (c) and (d). 

The complex behavior shown in the model HI2 im- 
ages is not seen in standard coronagraph images, and was 
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not found in the simula ted HI2 ima g es of a less struc- 
tured CME examined bv lLugaz et~all (|2005D . The reason 
the CME fades in and out is because of multiple cross- 
ings of the Thomson sphere by dense arcs that propagate 
thr ough the field of view of HI 2A as illustrated in Figure 
M (|Vourhdas k HowardI [2006h . Here, for the same time 
sequence as shown in Figure [51 we show the mass density 
in color on the equatorial (x — y) plane. We mark the inter- 
section of the Thomson sphere with the equatorial plane 
with a white circle. The position of STEREO A is marked 
along with black lines roughly showing the field of view 
of the HI2 coronagraph. Dense plasma is preferentially 
brightened in close proximity to the Thomson sphere. At 
t — 9.1 hours, a dense arc of plasma crosses the Thom- 
son sphere from the right hand side. This arc begins to 
move away from the Thomson sphere followed by a cavity 
at times t = 13.1 and 18.3 hours, which explains the dim- 
ming on the right hand side of the field of view. By time 
t = 22.3 hours, dense plasma on the Thomson sphere is 
found directly in front of the space craft on the left hand 
side, and far across the Thomson sphere on the right hand 
side. This two-sided distribution of plasma and its proxim- 
ity to the Thomson sphere coincides with the bright arcs 
seen in the left and right hand sides of the field of view 
seen in Panel (d) of Figure [H 

6. SUMMARY AND CONCLUSIONS 

We have preformed a detailed analysis of a 3D MHD 
time-dependent simulation of flux-rope-driven CMEs orig- 
inating from AR 10408 on October 26 and 28, 2003 
(|T6th et al. . 2007). The simulation was carried out with 
the SWMF in which propagation of the CMEs were fol- 
lowed from the low corona to beyond 1 AU. We have cre- 
ated synthetic coronagraph images of the October 28 event 
and made strict quantitative comparisons of these images 
with those made my LASCO C2 and C3. The simulation 
closely reproduces the appearance of the October 28 CME, 
both in the morphology and magnitude of Thomson scat- 
tered light through the fields of view of both LASCO C2 
and C3. Furthermore, the CME masses calculated from 
the synthetic coronagraph images and the LASCO obser- 
vations also are in close agreement. 

We have identified a shock front at the southern edge of 
the CME. A detailed analysis shows that the appearance 
of the CME-driven-shock in a synthetic coronagrpah im- 
age closely matches the same structure found the LASCO 
C3 image at 12:18 UT on October 28. Both the rigor of 
this quantitative comparison and the fidelity in reproduc- 
ing coronagraph images of CMEs is unprecedented. This 
is all the more remarkable in that there was no parame- 
ter fitting to produce the CME images, as only the CME 
speed at 20 Rs was used as a guide in specifying the fiux 
rope. 

The simulation allows the opportunity to view the Oc- 
tober 28 CME as it wo uld appear at larg e elongation 
in the HI2 coronagraph (|Lugaz et al.l 120051 ) . Here, we 
found complex time evolution of the CME in the HI2 
field of view as dense structures were temporarily high- 
lighted as they came in close proximity to the Thomson 
sphere. These results vividly illustrate the effects of CME 
pro pagation viewed at large elongation pointed out by 
l\^urlidas fc Howard (,2006,) . 



An aspect of the simulation that can be improved upon 
is the speed of the CME, which is a factor of two too fast 
inside of 5Rs, and does not recover the measured veloc- 
ity until 20 Rs. We conclude that most of the errors in 
CME speed originate from the ad hoc initiation with a 
fiux rope in a state of force imbalance. However, the ex- 
pansion of the flux rope, in the realistic 3D coronal model, 
is remarkable in its ability to reproduce the appearance of 
the CME in Thomson scattered white light. This result 
suggests that much of the density structure of the CME 
is determined by the plasma distribution in the corona 
surrounding the active region. 

Future models of CME initiation will be improved by 
the use vector magnetogram data rather than just fitting 
fiux ropes to the observed line-of-sight component of the 
field. For this event, vector magnetograms clearly show 
that the magnetic field of AR 10486 is very highly sheared 
with the field running nearly parallel to the neutr al line 
from which the CMEs erupt (|Liu &: Havashill2006l ). The 
evolution of magnetic shear is made manifest by observa- 
tions of photospheric proper motions, which show strong 
shear flows along t he magnetic neutr al line prior to CMEs 
(jYang et all 120041 : iDeng et al.ll2006f) . These observations 
strongly suggest that the magnetic shear is ultimately 
driving these CMEs. A physical explanation and simula- 
tion of such she ar flows was first provided respectively by 
[Manchester fc Low ( 2000) and Manchester (,2001) . which 
shows that the flows are driven by the Lorentz force that 
occurs as the magnetic field emerges in a stratified atmo- 
sphere. Furthermore, these Lorcntz-forcc-drivcn flows may 
persist until the magnetic field becomes so highly sheared 
that there is a loss of equilibrium p r oducing eruptions 
that ca n drive CMEs lManchesteil (|2003[) : [Manchester et all 
([2004c[ ). The need to resolve the photospheric pressure 
scale height has limited these simulations to Cartesian do- 
mains. Currently, we are adressing the challenge to include 
this self-consistent shearing mechanism in a global model 
of CME initiation that is guided by vector magnetic field 
observations. 
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